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We revisit the problem of the correlation gap in (n, n) armchair carbon nanotubes, that would be 
metallic in the absence of electron-electron correlations. We attack the problem in the context of a 
Hubbard model with on-site repulsion U only, and we show that the scaling of the gap as exp(— nf/f/) 
predicted by Balents and Fisher (Phys. Rev. B 55, R11973 (1997)), can only be valid if U is not 
too large, even for very large values of n. Using Hartree-Fock calculations and Renormalisation 
Group arguments we derive the scaling of the gap as a function of n for a given value of U . Possible 
applications for the magnitude of the correlation gap in armchair carbon nanotubes will be discussed. 

PACS numbers: 71.10.Pm, 71.20.Tx, 72.80.Rj 



I. INTRODUCTION 

In the last few years exciting experiments have been 
done on single wall carbon nanotubes (SWCNT). It has 
been shown that SWCNT can behave as a true quantum 
wire showing coherent transport [Q , Luttinger-liquid be- 
haviour was proposed |^ and there are indications for 
intrinsic superconductivity The typical dimensions 
of the CNT used in these experiments are a few nm for 
the diameter and several hundreds of nm for the length. 
Treating CNT in a single particle picture, the electronic 
properties depend strongly on chirality y. In particu- 
lar armchair CNT are metalic and their Fermi surface at 
half-filling consists of two points only. At this so-called 
Fermi points, bands with nearly linear dispersion rela- 
tions cross (Fig. 0). This property makes armchair CNT 
an ideal candidate for a one-dimensional quantum wire 
forming a Luttinger liquid 1^. 

In 1997 Balents and Fisher Q considered the problem 
of the correlation gap in half-filled (n, n) armchair carbon 
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FIG. 1: The energy bands near the Fermi level for armchair 
CNT at half-filling. 



nanotubes. By excluding all but the lowest bands and us- 
ing an on-site Hubbard interaction U , they mapped the 
problem onto a two-chain Hubbard model with an effec- 
tive interaction u„ = U jn. For this model it is known 
from renormalisation group (RG) calculations (see 
and references therein) that the functional dependence 
of the gap A on it„ and t is given by 



texp(— ct/M„) if M„ ^ i. 



(1) 
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In other words, (n, n) armchair CNT at half-filling with 
large enough n are predicted to be metallic for practical 
purposes since the correlation gap is exponentially small. 

In this paper we argue that the scaling law in Eqs. (|l|) 
can only be valid, even for large n, if at the same time 
the interaction strength I] is not too large. This can al- 
ready be seen from the ID limit, i.e. when n approaches 
infinity. In this limit we expect a metal-insulator transi- 
tion at some critical value \Jcr of the interaction strength 
[0, 1^ . Even at finite n the gap will be exponentially small 
only up to but will grow linearly for larger values of 
U . In the following we derive the scaling law of the gap 
as a function of n and U . This gap is determined from 
Hartree-Fock (H-F) calculations where we take into ac- 
count all the bands. We expect from the H-F method 
that it produces the correct functional dependence in the 
exponent of the correlation gap. This is the case for the 
Hubbard model in one dimension where the Bethe Ansatz 
solution is reproduced up to a prefactor On qualita- 
tive grounds we show how we can get this scaling law 
from a RG argument. 



II. THE CORRELATION GAP FROM H-F 
CALCULATIONS 

We are going to think about a SWCNT as a rectan- 
gular graphite monolayer rolled up into a cylinder. This 
is equivalent to considering a rectangular honeycomb lat- 
tice (Fig. H) with the appropriate periodic boundary con- 
ditions. All the SWCNT can be classified by their chi- 
rality vector Ch — nax + raa-2,, where a\ and 02 are 
the basis vector of the honeycomb lattice, while n and m 
are integers with to < n B. Ch determines into which 



2 



Ch 




-Tube axis 



FIG. 2: A rectangular honeycomb lattice of armchair type. 
The basis vectors are chosen to be ai = a/2(\/3, 1) and 
a2 = a/2{y/3,—l). a = 2.46A is the lattice constant of 
graphene. The chirality vector Ch points into the x direc- 
tion for armchair CNT. The tube axis is perpendicular to it 
in the y direction. 



direction the graphene layer is rolled up. As discussed 
in the introduction, the armchair CNT are of particular 
interest since at half-filling they are always metallic in 
the non-interacting electron approximation (see Fig. |^). 
Armchair CNT are characterized by chiral vectors of the 
form Ch = nai + na2, i.e. n — m. An example is shown 
in Fig. ^. In our calculations we take the tubes to be 
very long and we use also periodic boundary conditions 
at the ends of the tube. To determine the charge gap 
in (n, n) armchair nanotubes, we consider the following 
model with nearest neighbor hopping and on-site Hub- 
bard interaction at half-filling: 



H = 



(2) 



a is the spin index and i , j sum over the sites of a rect- 
angular armchair-type honeycomb lattice with periodic 
boundary conditions, cj^. are the fermion creation oper- 
ators and riia = cl^Ci^- The hopping integrals tij are re- 
stricted to nearest neighbors and are taken to be equal to 
a single t for all the hoppings. In the case of half-filling 
and on bipartite lattices with non-frustrated hopping, 
the H-F Hamiltonian is ||] 



H 



h.c. 



E/ ^i3^\<y^ 

+U ^ + {n^^)n^^ - {ni^){n^l)] (3) 



with the expectation values given by 

(rw) = 1/2 (l + m(-l)U|,) 



(4) 



where X^a = 1, ^a-a = — 1 and m = \ {ni-^) — {ni\)\. This 
Hamiltonian can be diagonalized by a Bogoliubov trans- 
formation and it yields the self-consistent H-F equation 
for the sublattice magnetisation: 
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FIG. 3: Numerical results for the correlation gap: (a) The 
gap for a (5, 5) tube as a function of U . The solid line is the 
H-F result. The dashed line is the fit to cU/texp{—Xt/U) (c, 
A are constants), (b) The gap as a function of the CNT size 
n at U — 2t. The H-F results (o) are fitted to cexp(— An) 
(sohd line) (c) The same plot for U = (/i^^ = 2.23t. The data 
fits the function c/n. 



N is the total number of sites and e(fc) is the tight binding 
dispersion relation for a single graphite layer given by: 

e2(fc) ^ + 2t'^ [cos{k ■ ai) + cos{k ■ a2) 

+ cos{k ■ (ai - 02))] (6) 

Finally the gap is obtained from the sublattice magneti- 
sation: 



A = min ^/lPm^ + 4e^k) = Um (7) 

k 

In Fig. ^ the numerical results for the correlation gap 
are shown. We have two parameters to vary, the tube 
diameter which is proportional to n and the interac- 
tion strength U. We tried to fit the H-F results to 
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cU/texp{—Xt/U) when we varied U and to cexp(— An) 
for a variation of n. We see that these fits break down at 
a critical value of the interaction strength = 2.23t. 
There is no exponential decay beyond C/^^. In particu- 
lar for U — we can fit the data to c/n, a power law. 
All our calculations can be summarized in the following 
scaling law 

A/t = l/nexp{-an{t/U -t/U^/)} (8) 

with a = 5.44. is identical to the critical value 

in the H-F approximation to open a gap in the two- 
dimensional honeycomb lattice. We see that if U is ap- 
proaching the exponent is going to vanish and we 
are left with a power law A = t/n. 



III. RG~ARGUMENT 
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FIG. 4: The correlation gap as a function of interaction 
strength for (5,5) CNT (dot-dashed line), (10,10) (dashed 
line) and the 2D-limit (solid line). 



To compute the charge gap from a RG-calculation we 
need to compute the RG-equation up to third order in 
the running coupling constant jl^. For armchair CNT 
the effective coupling constant is given hy g = U/n ||] 
and we denote by g the renormalized coupling constant 
corresponding to the momentum cutoff D. The third 
order equation is [ndl 



dlnL> 
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with Vc = vp + U/{2Trn). The RG flow has to be stopped 
at the strong coupling region. This should occur when 
the momentum cutoff D corresponds to the energy scale 
of the charge gap Ac ||l^ . The flow of the running cou- 
pling constant has to be stopped at the first relevant 
energy scale. In the usual one-band picture this energy 
scale is just the kinetic energy ~ t. But in the present 
problem there is another energy scale entering before: It 
is the gap to the next band which is of order t/n . This 
argument leads to g(Ac) ~ t/n. Now we can integrate 
out the degrees of freedom far from the Fermi level, and 
we get 



d(lnl)) 



9(Ae) 



dg { (l ^ 



(10) 

where Dq is the initial cutoff with g{Do) = g ^ U /n. To 
leading order in g and g{Ac) we obtain 
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(11) 



Identifying the relevant energy scales we can write the 
charge gap as 

Ac -exp(-/3t/([//n) -f 7n). (12) 
We note that this argument leads to the H-F result. 



IV. DISCUSSION 

Fig. 1^ summarizes our results obtained from H-F cal- 
culations and the RG argument on the correlation gap in 
armchair CNT. We observed that there is an exponen- 
tially small gap only if [/ <C , even for large n. 
is a critical value for the interaction strength. For values 
U < U^^ the correlation gap follows an exponential scal- 
ing law. For U > the gap is growing linearly. The 
actual value of the hopping integral t in SWCNT is about 
2.4 eV |l^. This value is large and as a consequence we 
observe in the case U < U^^ a finite gap which is of 
order meV. As an example we look at a (10, 10) CNT 
assuming an interaction strength of ?7 = 2t. In this case 
we obtain for the gap a value of 0.006i which corresponds 
to an energy of 14 meV and a temperature of 160 K. For 
smaller CNT the gap is even larger. 
These results are based on two approximations, namely 
the description of electron-electron interactions in terms 
of a Hubbard model with only on-site repulsion, and the 
calculation of the gap within Hartree-Fock theory. Let 
us now critically review these approximations. 
Hubbard model: It is by now well accepted that metallic 
nanotubes behave as Luttinger liquids with an exponent 
controlled by the long-range tail of the Coulomb repul- 
sion and by the finite length of the tube However, 
when it comes to calculating the correlation gap in half- 
filled systems, the short-range part of the correlations 
plays a dominant role. This is most easily seen in the 
atomic limit where the hopping integrals are assumed to 
be vanishingly small. So let us assume for a moment that 
the system is described by the Hamiltonian: 



H 



U ^n^^n,i 



^V^jn^nj 



(13) 



where = n^i -I- ^ is the total density, and where Vij 
is the long-range part of the repulsion that does not need 
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to be specified. According to Mott's prescription to eval- 
uate the gap one has to compare the energy of a 
uniform configuration with that of a configuration with 
a doubly occupied site and a hole far apart from each 
other. The only difference comes from the energy of the 
electron that has been moved, and the energy increase 
is precisely equal to U since it still interacts with fur- 
ther neighbours in the same way. So, even in the pres- 
ence of long-range Coulomb repulsion, the value of the 
charge gap is controlled by the on-site repulsion U in the 
atomic limit. Of course, away from the atomic limit, the 
long-range part of the Coulomb interaction will play a 
role. To get a quantitative estimate of the charge gap 
in that case is a very difficult problem though which has 
not been solved even in the simplest case of a pure one- 
dimensional model, but if anything the longer range part 
of the Coulomb repulsion is expected to reinforce the ten- 
dency to localize the charge, hence to increase the charge 
gap. So to use a simple Hubbard model with only on-site 
repulsion to calculate the correlation gap of a half-filled 
system is a reasonable assumption, and the value is prob- 
ably an underestimate of the actual gap in the presence 
of the long-range part of the Coulomb repulsion. For 
the present purpose this is all we need since our main 
conclusion is to argue that the gap might be larger than 
previously assumed. 

Hartree-Fock approximation: In purely one-dimensional 
systems, Hartree-Fock is known to reproduce correctly 
the exponential form of the gap as a function of U , but 
the prefactor of the exponential is wrong, and the gap is 
overestimated by a factor \JTJ\J . This is only a problem 
for very small values of U /t, but in the range of interest to 
us, namely not too far from Ucn this is not a problem any 
more. In two dimensions, very little is known precisely re- 
garding the value of the correlation gap for models with- 
out perfect nesting. Luckily enough, the Hubbard model 
on the honeycomb lattice has been extensively studied by 



Monte Carlo simulations, and although the gap could not 
be calculated, the critical value Ucr has been determined 
by looking at the development of magnetic long-range or- 
der 10, ^, and the value is around Ucr — 3.6t. This is 
clearly larger than the Hartree-Fock result = 2.23t, 
but the order of magnitude is the same. Since the gap 
is anyway expected to grow linearly with U beyond the 
critical value, the main result of Hartree-Fock that the 
gap in the nanotube geometry could take large values if 
U is not too far below Ucr is expected to remain true 
beyond Hartree-Fock. 

The relevance of this calculation depends crucially on 
the actual value of U in nanotubes. What we need here 
is the atomic value for carbon. The fact that graphite, 
a system of very weakly coupled honeycomb planes, is 
a semi-metal and not an insulator just tells us that U 
does not exceed Ucr- Besides, values often quoted for 
fuUerenes are of no help since they concern the molecule 
Ceo and not the atomic value for carbon. The most rel- 
evant value which we could find in the literature was 
for polyacetylene where an on-site repulsion of 5-10 eV 
is generally used Q. So we believe that values in the 
range of 2t to 4t are to be expected. In particular, a value 
of U near the critical one seems plausible. Thus if in real 
CNT U < Ucr, a correlation gap of a few meV has to be 
present and a gap of such magnitude could in principle 
be observed experimentally. 

This of course will only be true for SWCNT at half- 
filling. Currently available CNT seem to have some self- 
doping which shifts the Fermi-level to an energy where 
several bands are cut. This turns the CNT metallic in- 
dependently of their chirality. There are ongoing experi- 
mental efforts however identifying the nature of the dop- 
ing and producing undoped samples, and we are confident 
that it will be soon be possible to test the conclusions of 
the present work. 
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